Binary upscaling—the role of connectivity and a new formula     

Christen Knudby  Jesús Carrera  Jim D. Bumgardner  Graham E. Fogg 《Advances in water resources》2006

Although recognized as important, measures of connectivity (i.e. the existence of high-conductivity paths that increase flow and allow for early solute arrival) have not yet been incorporated into methods for upscaling hydraulic conductivities of porous media. We present and evaluate a binary upscaling formula that utilizes connectivity information. The upscaled hydraulic conductivity (K) of binary media is determined as a function of the proportions and conductivities of the two materials, the geometry of the inclusions, and the mean distance between them. The use of a phase interchange theorem renders the formula equally applicable to two-dimensional media with inclusions of low K and high K as compared with the matrix. The new upscaling formula is tested on two-dimensional binary random fields spanning a broad range of spatial correlation structures and conductivity contrasts. The computed effective conductivities are compared to what is obtained using self-consistent effective medium theory, the coated ellipsoids approximation, and to a streamline approach. It is shown that, although simple, the proposed formula performs better than available methods for binary upscaling. The use of connectivity information leads to significantly improved behavior close to the percolation threshold. The proposed upscaling formula depends exclusively on parameters that are obtainable from field investigations.  相似文献   

Middle Miocene to Quaternary primary basalt and high magnesian andesite magmas of North Hokkaido, Japan: Source mantle characteristics and degrees of partial melting   总被引：1，自引：0，他引：1  

Hiroyuki  Ishimoto  Kenji  Shuto  Yoshihiko  Goto 《Island Arc》2006,15(2):251-268

Abstract   Middle Miocene to Quaternary primitive basalts and high magnesian andesite (HMA) in North Hokkaido resulted from three periods of intense volcanism; early-stage (12–10 Ma), middle-stage (9–7 Ma) and late-stage (3–0 Ma). Based on the chemical compositions of olivines and chromian spinels and bulk chemistry of the primitive rocks, we examined depths of segregation of the calculated primary magmas and the degrees of partial melting of the source mantle. In the context of asthenospheric mantle upwelling, petrological data from the present study can be accounted for by the secular change in the depth of magma segregation from the upwelled asthenospheric mantle, which is composed of fertile peridotite. Thus, the early-stage primary magmas were generated by higher degrees of partial melting of the shallower part of hot asthenospheric mantle, whereas the middle- and late-stage primary magmas resulted from lower degrees of partial melting of a deeper part of the asthenospheric mantle. The early-stage HMA magma was generated by partial melting of the remnant subcontinental lithospheric mantle composed of refractory peridotite. This melting might have resulted from an increased geothermal gradient caused by upwelling of hot asthenosphere.  相似文献   

Rheological and kinematical responses to flow of two-phase rocks     

Yoshi-Taka Takeda  Albert Griera 《Tectonophysics》2006,427(1-4):95

Numerical simulations have been performed to investigate the strain-dependent behaviour of rheological and kinematical responses to flow of two-phase rocks using the commercial finite-difference program FLAC2D. It was assumed that the two phases have Maxwell rheology. Plane strain and velocity boundary condition, which produces a simple shear deformation, were also assumed. Two types of geometries were considered: strong phase supported (SPS) and weak phase supported (WPS). We calculated strain-dependent variations of effective viscosity and partitioning of strain rate, vorticity and kinematic vorticity number during deformation in both SPS and WPS structure models.The results show that the strain-dependent behaviour is largely influenced by the geometry of the composite. SPS models show both strain hardening and strain softening during the simulations, with strain hardening preceding strain softening. A critical shear strain is necessary to begin the strain softening behaviour. Strain hardening and strain softening are accompanied by a reduction and an increase of the partition of strain rate into the weak phase, respectively. On the other hand, WPS models show only weak strain hardening and strain softening, being the strain-dependent behaviour close to a steady state flow. In addition, the following results are obtained on vorticity and kinematic vorticity number; (1) in both SPS and WPS models the partition of vorticity into weak phase increases with progressive shear strain, i.e. the strong phase becomes less rotational, (2) in SPS models weak inclusions changes from sub-simple shear to super-simple shear with progressive strain, whereas the strong matrix changes from super-simple shear to sub-simple shear, (3) in WPS models the strong inclusions with high viscosity contrasts are less rotational but can be in super-simple shear condition to high strains.The observed strain-dependent behaviours have been compared with previous proposed analytical models. The degree of agreement is variable. Balshin and Ryshkewitch–Duckworth models are only applicable to SPS models. Ji-generalized mixture rule model is applicable to both models.The results suggest that polyphase rocks with SPS structure during ductile shear deformation respond as strain softening materials, after an initial strain hardening stage that may drive to the strain localization into the material.  相似文献   

A k - *epsiv Turbulence Closure Model For The Atmospheric Boundary Layer Including Urban Canopy     

Vu Thanh Ca  Yasunobu Ashie  Takashi Asaeda 《Boundary-Layer Meteorology》2002,102(3):459-490

A numerical model for the computation of the wind field,air temperature and humidity in the atmospheric boundary layer (ABL) including the urbancanopy was developed for urban climate simulation. The governing equations of the modelare derived by applying ensemble and spatial averages to the Navier–Stokes equation, continuityequation and equations for heat and water vapour transfer in the air. With the spatial averagingprocedure, effects of buildings and other urban structures in the urban canopy can be accounted for byintroducing an effective volume function, defined as the ratio between the volume of air in acomputational mesh over the total volume of the mesh. The improved k - model accounts for the anisotropyof the turbulence field under density stratification. In the improved k - model, the transportof momentum and heat in the vertical direction under density stratification is evaluated based onthe assumption of a near-equilibrium shear flow where transport effects on the stresses andheat fluxes are negligible. The heating processes at surfaces of buildings and ground are alsomodelled. The comparison of the computational results obtained with the present modeland existing observational data and numerical models shows that the present model is capableof predicting the structure of turbulence in the urban canopy layer under density stratification.Numerical experiments with the new model show that the flow behaviour of the air in the urbancanopy layer is strongly affected by the existence of buildings and density stratification.  相似文献   

The Relationships between Effective Porosity,Uniaxial Compressive Strength and Sonic Velocity of Intact Borrowdale Volcanic Group Core Samples from Sellafield     

D.?C.?EntwisleEmail author  P.?R.?N.?Hobbs  L.?D.?Jones  D.?Gunn  M.?G.?Raines 《Geotechnical and Geological Engineering》2005,23(6):793-809

The effective porosity, saturated sonic velocity and saturated uniaxial compressive strength were determined on a large number of Borrowdale Volcanic Group volcaniclastic core samples from three boreholes at Sellafield, Cumbria. The work formed part of the UK Nirex Limited site investigation into whether the Sellafield area could be suitable as a repository for intermediate and low level radioactive waste. Most of the intact samples were of low to very low effective porosity, had a high sonic velocity and were very strong to extremely strong. However, a proportion of values deviated significantly from this. Bivariate analysis showed a negative relationship exists between sonic velocity and effective porosity. The cross plots of these two parameters with uniaxial compressive strength showed a wide range of strength values for samples of low effective porosity and high sonic velocity. Six failure types were identified during the uniaxial compressive strength tests. The strongest samples tended to fail through the matrix and the weakest rock samples tended to fail through haematized material or along haematized veins. Effective porosity and sonic velocity measurements could not distinguish between those samples that failed through the matrix and those that failed along discrete narrow veins. The presence of narrow haematized veins has a major effect on the intact rock strength.  相似文献   

不同碳源强化地下水中生物脱氮模拟试验研究   总被引：12，自引：0，他引：12  

张云  张胜  杨振京  张凤娥 《地理与地理信息科学》2003,19(1):66-69

地下水中NO3^-N污染是普遍存在的环境问题,生物脱氮作用是去除该污染的主要机制,而生物脱氮菌群和营养碳源又是这种作用进行的主要限制性因子,该文应用这一理论将人工接种驯化并优势培养制备的生物脱氮菌剂,分别与不同种类不同配入量比的营养碳源物质一起施用,进行系列污染水体的生物脱氮模拟试验研究,以确定用于治理地下水中NO3－N污染的微生物菌剂和促进生物脱氮作用的营养碳源种类及其最佳配入量比值,探索修复治理地下水中大面积NO3－N污染的方法。  相似文献   

The underground application of Magnetic Resonance Soundings   总被引：4，自引：0，他引：4  

J.M. Greben  R. Meyer  Z. Kimmie 《Journal of Applied Geophysics》2011,75(2):220-226

The potential application of MRS technology in locating waterbearing fractures in underground mines is studied. The determination of the presence of water ahead of mining is important to prevent accidents and to ensure higher efficiency in mining operations. In the usual surface based measurements, with horizontal loop and water layer, the geometry of the problem can be summarized by the value of the inclination of the Earth magnetic field. For MRS measurements under the geometric conditions associated with underground mining, where the loop is non-horizontal, the geometry can be described in an effective inclination that can be expressed in terms of the Earth magnetic inclination and declination, together with two further parameters that characterize the orientation of the mine wall. We examine the consequences of the different geometries on the MRS signal. Since the loop size is severely restricted in underground conditions, the feasible target depth is also severely limited. The consequences of the fractured hard rock aquifer conditions, typical of deep mining or tunneling environments, are also examined. The overall conclusion is that in underground MRS applications the signal strength is too small to enable the practical identification of fractures containing large volumes of water ahead of the mining face.  相似文献   

Truncated multiGaussian fields and effective conductance of binary media     

Sean A. McKenna  Jaideep RayYoussef Marzouk  Bart van Bloemen Waanders 《Advances in water resources》2011,34(5):617-626

Truncated Gaussian fields provide a flexible model for defining binary media with dispersed (as opposed to layered) inclusions. General properties of excursion sets on these truncated fields are coupled with a distance-based upscaling algorithm and approximations of point process theory to develop an estimation approach for effective conductivity in two-dimensions. Estimation of effective conductivity is derived directly from knowledge of the kernel size used to create the multiGaussian field, defined as the full-width at half maximum (FWHM), the truncation threshold and conductance values of the two modes. Therefore, instantiation of the multiGaussian field is not necessary for estimation of the effective conductance. The critical component of the effective medium approximation developed here is the mean distance between high conductivity inclusions. This mean distance is characterized as a function of the FWHM, the truncation threshold and the ratio of the two modal conductivities. Sensitivity of the resulting effective conductivity to this mean distance is examined for two levels of contrast in the modal conductances and different FWHM sizes. Results demonstrate that the FWHM is a robust measure of mean travel distance in the background medium. The resulting effective conductivities are accurate when compared to numerical results and results obtained from effective media theory, distance-based upscaling and numerical simulation.  相似文献   

某核电站导流堤地基液化及地震永久变形分析     

赵杰  郑洋  王桂萱  周金领 《地震工程学报》2011,33(Z1)

在地震反应分析中常采用总应力分析法,但总应力法没有考虑孔隙水压力的变化规律和液化随时间的发展过程。基于二维有效应力动力分析方法,结合Biot动力固结理论,采用自行开发的有效应力动力分析程序对某核电站导流堤地基进行液化分析,给出在地震作用下砂土层的液化范围,并计算出永久变形。所得结论可以给类似工程提供理论指导。  相似文献   

Behavior of vortex-induced vibration of a circular cylinder near a deformable wall with two degrees of freedom in steady flow     

杨兵  马建林  崔金声  徐万海 《中国海洋工程》2011,25(4):549-564

The behavior of vortex-induced vibration of a two-degree-of-freedom cylinder near a deformable wall in steady flow is investigated experimentally.The typical phenomenon of the two-degree-of-freedom cylinder’s VIV is discussed.The influences of initial gap between the cylinder and the wall on the dynamic responses of the cylinder are analyzed.The comparison is made about dynamic responses of the cylinder with one and two degrees of freedom.Experimental results show that the vibration of the cylinder near a deformable wall with a small value of initial gap-to-diameter ratios can generally be divided into two phases.The initial gap-to-diameter ratios have a noticeable influence on the occurrence of transverse vibration.The transverse maximum amplitude of the cylinder with two degrees of freedom is larger than that of the cylinder with one degree of freedom under the condition with the same values of other parameters.However,the vibration frequency of the cylinder for the two degrees of freedom case is smaller than that for the one degree of freedom case at the same value of Vr number.  相似文献